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ABSTRACT 

Black hole masses for samples of active galactic nuclei (AGN) are currently es- 
timated from single-epoch optical spectra. In particular, the size of the broad-line 
emitting region needed to compute the black hole mass is derived from the optical 
or ultraviolet continuum luminosity. Here we consider the relationship between the 
broad-line region size, R, and the near-infrared (near-IR) AGN continuum luminosity, 
L, as the near-IR continuum suffers less dust extinction than at shorter wavelengths 
and the prospects for separating the AGN continuum from host-galaxy starlight are 
better in the near-IR than in the optical. For a relationship of the form R oc L", we 
obtain for a sample of 14 reverberation-mapped AGN a best-fit slope of a = 0.5 ±0.1, 
which is consistent with the slope of the relationship in the optical band and with the 
value of 0.5 naively expected from photoionisation theory. Black hole masses can then 
be estimated from the near-IR virial product, which is calculated using the strong and 
unblended Paschen broad emission lines (Paa or Pa/3). 
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1 INTRODUCTION 

The discovery of tight correlations between a galaxy's 
central black hole mass and the luminosit y and veloc- 
ity dispersion of its stellar bulge jM agorria n et al. 1 1 19981 : 
iGebhard t et al.' 2OO0I: iFerrarese fc Merrittll200(il ') is expected 
to strongly constrain how galaxies form and grow over cos- 
mic time. Therefore, much effort goes in par ticular into mea- 
suring the rate of black hole growth (e.g., Yu fc Tremaind 



l2002l : THeckman et al.|[20o3 : iKellv et alj|2009l ). Since this re- 
quires both large samples of galaxies with easily obtainable 
black hole mass estimates and sources that probe the highest 
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redshifts, such studies rely heavily on active galactic nuclei 
(AGN). 

In AGN, black hole masses can be readily estimated 
from single-epoch optical spectra. Assuming that the dy- 
namics of the broad-emission line gas seen in these sources 
is dominated by the gravitational force of the black hole, 
one can use the virial theorem to calculate black hole masses 
from only two observables, the velocity dispersion and ra- 
dial dist ance of the emitt ing gas from the central source 
(e.g., Peterson et al.l [20041 ') . The velocity dispersion can be 
obtained from the widths of the Balmer hydrogen lines 
(H/3 or Ha) and the continuum luminosity at (rest-frame) 
5100 A serves as a surrogate for the size of the broad-line 
regi on (BLR ) . The l atter finds its j ustific ation in the work 
of PetersonI (|l993l ). IWandel et all (|l999l ) and iKaspi et all 
(|2000l ). who showed that emission-line lags, and so broad- 
emission line radii, derived from optical reverberation map- 
ping studies correlate with the optical luminosity (of the 
ionising component) largely as expected from simple pho- 
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Table 1. Reverberation-Mapped Sample of AGN 



Object Name 
(1) 


Rap 
[It-days] 
(2) 


log I^Llpm 

[erg s~i] 
(3) 


log ^L^^^ 
[erg s 1] 

(4) 




307I91 


45.81±0.Ui 


45.d4±U.1U 


Mrk876 


40±15 


44.77 


44.55 


PG0844+349 


32+" 

— io 


44.28±0.02 


44.28±0.10 


MrkllO 




43.69±0.08 


43.66±0.05 


Mrk509 


80+^ 


44.05±0.03 


44.07±0.04 


Arkl20 




44.26±0.06 


43.70±0.10 


Mrk817 


22+_l 


43.71±0.03 


43.64±0.01 


Mrk290 


9±1 


43.43±0.04 


43.22±0.12 


Mrk335 


16l^ 


43.49±0.05 


43.35±0.07 


Mrk79 




43.30±0.04 


43.20±0.05 


NGC5548 


18.0 ±0.6 


43.12±0.05 


42.97±0.09 


NGC7469 


4 5+0.7 


43.45 


43.10 


NGC4593 


3.7±0.8 


42.73±0.02 


42.50±0.06 


NGC4151 


7±1 


42.13±0.06 


41.84±0.05 



The columns are: (1) object name; (2) radius of the 
H/3 broad-emission line region (in light-days) (Ha for 
PG 0844-1-349); average integrated 1 /xm continuum lu- 
minosity (4) as observed and (5) fiux-calibrated relative 
to best- weather epoch and host galaxy-subtracted. Errors 
represent la uncertainties. For sources writh multiple radius 
measurements we list the weighted average (in boldface). 



toionisation arguments. This correlation is now referred to 
as the radius-luminosity (R-L) relationship. 

We consider here the application of this technique in 
the near- infrared (near-IR), which potentially aflfords some 
advantages over the optical (and ultraviolet). First, the ob- 
served optical continuum can be severely contaminated by 
host galaxy starlight if a large slit is used (as is often the case 
in reverberation studies), especially in low-redshift sources 
that have a weak AGN or sources with a luminous stellar 
bulge, and so does not give directly the ionising flux. High- 
resolution, deep optical images are then required to estimate 
the host galaxy starlight encl osed in the aperture in order 
to correct the optical spectra (|Bentz et al ] |2006l . l2009l ). Sec- 
ondly, the optical hydrogen broad-emission lines, in partic- 
ular H/3, are strongly blended with other species, making 
it difficult to reliably measure the line width. Thirdly, all 
optical measures can sufTer from dust extinction. 

Here we present the radius-luminosity relationship in 
the near-IR, which, in conjunction with the width of the 
strong and unblended Paschen broad emission lines Paa and 
Pa/3, can be used to estimate AGN black hole masses. In 
Section 2, we introduce the data and present the near-IR 
R-L relationship. Its application is discussed in Section 3. 
In Section 4, we present our conclusions. Throughout this 
paper we have assumed cosmological parameters Ho — 70 
km s"^ Mpc"^ Qm = 0.3, and = 0.7. 



2 THE NEAR-INFRARED R-L RELATIONSHIP 

We have recently shown l|Landt et al.ll201lh that in broad- 
emission line AGN the accretion disc spectrum, which is 
believed to be the main source of ionising radiation, extends 
well into the near-IR and still dominates the continuum at 



100 



10 




log (i^LiMm/lO"') ferg/sl log (i^Lt™/lO«) ferg/s] 



Figure 1. The radius of the H/3 broad-emission line region (in 
light-days) versus the average 1 fim continuum luminosity as ob- 
served (left panel) and flux-calibrated relative to best-weather 
epoch and host galaxy-subtracted (right panel). For sources with 
multiple radius measurements we plot the weighted average. The 
dashed, solid and dotted lines indicate the correlations obtained 
with the BCES, FITEXY and GaussFit routines, respectively. 



~ 1 /im. Therefore, a single-epoch near-IR spectrum can in 
principle be used to estimate the broad-emission line radius. 

In Fig. [1] we verify this conjecture by plotting the ra- 
dius of the H/3 broad-emission line region (Rnp) versus the 
integrated 1 fim continuum luminosity for the reverberation- 
mapped AGN in our sample (14 objects; see Table [1]). Val- 
ues for the radius of the H/3 broad-line re gion are from 
iBentz et al] (|2009l ) and for Mrk 290 from iDennev et all 
( 20101 ). The near-IR measures are based on data obtained 
with the SpeX spectrograph (|Ravner et all |2003| ) at the 
NASA Infrared Telescope Facility (IRTF) on Mauna Kea, 
Hawai'i. The excellent atmospheric seeing at this observ- 
ing site allowed us to use a relatively narrow slit of 0.8", 
which excluded most of the host galaxy starlight. The only 
exceptions were Mrk 590 and NGC 3227, two reverberation- 
mapped AGN in our sample that were found to be in a very 
low state with the continuum s trongly dom i nated by host 
galaxy emission (see Fig. 6 in iLandt et al.l I2OIII ) . There- 
fore, they are omitted from the present study. We observed 
sources on average twi ce within a period of ~ 3 years 
(|Landt et al.ll2008[ I2OIII ). Here we use the mean integrated 
1 continuum luminosity and the error on the mean. We 
consider both the observed values (vLif^^) and the pure 
AGN values (lyLt^^). The latter were derived by perform- 
ing a relative flux-calibration using the [S III] A9531 nar- 
row emission line in the best-weather epoch and correcting 
for host-galaxy s tarlight using the Hubbl e Space Telescope 
(HST) images of lBentz et ai] |200d [2OO9I ) and for Mrk 290 
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a recent, unpublished HST im age (Bentz e t al. i n prep.). 
Further details are provided bv lLandt et al.l (|201ll ). 

We have performed linear fits of the form 
loR Rub ~ K + a \o^(u L^ / 10^ ^ ) following the approach 
of iBentz et al.1 l|2006l ) and iBentz et al.l l|2009l ). In par- 
ticular, we have used the three fitti ng routines BCES 
l|Akritas fc Bersha dv 1996 ), FITEX Y (Press et al.l |2007| ) 
and GaussFit (^McAr thur et al ][l99i) that can incorporate 
errors on both variables and, except for GaussFit, allow us 
to account for intrinsic scatter. Note that accounting for 
intrinsic scatter has the effect of increasing the weight given 
to data points with the largest errors, which is preferred 
if the i ntrinsic dispersion is larger than the measurement 
errors (ITremaine et al. I l2002i l. Roughly half of our sample 
(6/14 sources) has multiple measurements of 7?h/3 and, 
therefore, we have considered the following two cases: (i) 
using the average derived from all measurements for a 
particular source and weighted by the mean of the positive 
and negative errors (i.e., weighted averages) and (ii) using 
Monte-Carlo (MC) techniques to randomly sample _Rh/3 
from the individual values for each object. Two sources 
have only one epoch of near-IR data and in these cases we 
have assigned to the logarithmic luminosity the average 
error of the sample. We obtain significant correlations in all 
cases. Table[2]lists the results, which are shown for the case 
of the weighted averages in Fig. [T] Note that the near-IR 
continuum luminosities and the optical broad-emission line 
radii are not measured simultaneously, which is expected to 
increase the scatter in their relationship. 

All six cases for both total and AGN 1 ^m coritirmum 
luminosity give similar values for the slope and the intercept. 
The slope is in most cases ~ 0.5 ±0.1 and always consistent 
with a value of 0.5. Simple photoionisation arguments sug- 
gest that a given emission line will be produced at the same 
ionising flux in all AGN and therefore ii cx L^^^. Our results 
are si milar to the best-fit slope of 0.521q q7 found in the opti- 
cal bv lBentz et al.l l|200S l. who used host galaxy-subtracted 
fluxes of 34 reverberation-mapped AGN (more than twice 
our sample size). 



3 DISCUSSION 

Given the importance of the R-L relationship for AGN black 
hole mass determinations, alternatives to the optical con- 
tinuum luminosity are already being explored, e.g., the X- 
ray l u minosity and broad H/9 fine luminosity (|Kaspi et al.l 
l2005l : iGreene et allbOlOl). the [O il l] A5007 A emission- 

t he Pa g and Pa/3 
and 



line luminosity I Greene et aL 20101) 



emission-line luminosities I Kim et al] |2010| ). and the 
[Ne v] A 14.32 /xm and JO iv ] A25.89 nm emissio n-line lumi- 
nosities ()Dasvra et al.ll2008l : [Greene et al.ll2010l) . The 1 fim 
continuum luminosity is an efficient alternative that con- 
nects directly to the spectrum of the ionising source. 

Near-IR spectroscopy is now available at excellent ob- 
serving sites that regularly achieve subarcsecond seeing. 
This means that the host galaxy flux contribution in the 
near-IR can be reduced to a minimum. For example, a slit 
width (and seeing) of 0.8" in the near-IR relative to one of 
3", which is typical of most optical telescope sites, leads to 
a reduction of the extraction aperture (and so of the host 
galaxy flux contribution, assuming constant surface bright- 




Figure 2. The aver age integrated 5100 A continuum luminos- 
ity from lBentz et al.l (2009) versus the average integrated 1 fim 
continuum luminosity. Open and filled circles indicate total and 
host galaxy-corrected optical values, respectively. The host galaxy 
correction to the optical luminosities has the effect of changing 
the observed correlation slope from 0.91 ± 0.04 (dashed line) to 
1.03 ±0.07 (solid line). 



ness) by a factor of ~ 14. The case is exacerbated if we 
consider the much larger apertures employed in optical re- 
verberation studies (5 — 10"). Additionally, the seeing is im- 
proved at longer wavelengths, thus reducing the extraction 
aperture further (e.g., by a factor of ~ 1.3 at 1 /im rela- 
tive to 5100 A). Adaptive optics should make it possible 
to realise even greater gains in the near-IR relative to the 
ground-based optical. 

As our results show, the correlations for the total and 
AGN 1 fj,m continuum luminosity are consistent with each 
other and so independent of the host galaxy flux subtrac- 
tion, although the slope is more often ~ 0.5 and the scatter 
slightly smaller for the latter. Performing the same flts for 
our sample (14 sources) using the optical data of lBentz et al.l 
(2009. ) with the host galaxy flux subtracted gives results 
that are more similar to the total rather than to the AGN 
near-IR spectra (see Table [2|. However, if we consider their 
entire optical sample (34 sources), we obtain results similar 
to the AGN near-IR case. In this respect, we note that the 
intercepts of the optical and near-IR R-L relationships are 
similar, a result which is not unexpected if the two AGN 
luminosities sample the canonical accretion disc spectrum 



.1/3 



). In this case, the difference in integrated lu- 



minosity would be ~ 0.4 dex, well within the scatter of 
~ 0.5 dex of the near-IR and optical relationships for our 
sample (see Table[2|. Fig.[2]further illustrates the efficacy of 
our approach to minimising the host-galaxy contamination 
of the near-IR AGN continuum. While the near-IR lumi- 
nosities we measure correlate well with the measured optical 
luminosities (open circles), with a slope of 0.91 ± 0.04, they 
correlate even better once the host-galax;y contributions to 
the optical luminosities have been removed. The slope of the 
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Table 2. Best-Fits for the Relation logi?H/3 = K + o log(i^L^/10'*3) 





vL 


Ipm (14 obj.) 




uL 


(14 Obj.) 




uL^'^^. (14 obj.) 
5100A ^ ' 




uL^'^^, (34 obj.) 

5100A ^ ' 




Type 


K 


a 


S* 


K 


a 


S* 


K 


a 


S* 


K 


a 


S* 


BCES 


Avg 
MC 


0.91±0.11 
0.90±0.13 


0.56±0.10 




1.02±0.09 

1 01+"" 


0.55±0.10 




0.91±0.12 
0.91±0.14 


0.56±0.12 
0.56±0.13 




1.01±0.07 
0.98±0.10 


0.52±0.05 
0.52±0.07 




FITEXY 


Avg 
MC 


0.82±0.10 
0.81±0.11 


0.52±0.10 
0.53±0.11 


53 
52±2 


0.95±0.08 
0.93±0.09 


0.52±0.08 
0.52±0.09 


48 
48±2 


0.84±0.09 
0.82±0.11 


0.53±0.10 
0.53lO;l? 


50 
51±2 


0.91±0.06 
90+° °* 


0.54±0.05 
0.54±0.06 


40 
41±1 


GaussFit 


Avg 
MC 


l.OOitO.lO 
0.95±0.14 


0.54±0.10 
0.56±0.12 




1.09±0.08 
1.07±0.11 


0.54±0.09 

54+0'09 




0.97±0.09 
0.92±0.12 


0.56±0.09 
0.58±0.10 




1.01±0.05 
0.96±0.08 


0.53±0.04 
0.56±0.05 





* Scatter calculated as the percentage of the log Rhi3 value that, when added in quadrature to the error value, gives x2 = 1 



relationship between the optical AGN luminosities and the 
near-IR luminosities is 1.03 ± 0.07. Clearly the narrow-slit 
near-IR luminosities are dominated by the AGN component. 

As we noted earlier, attenuation by dust is much re- 
duced in the near-IR compared to the optical, so the 1 pm 
continuum luminosity may be particularly useful as an al- 
ternative for dust-obscured AGN. A modest intrinsic ob- 
scuration of a hydrogen column density of A^h ~ 10^^ cm"'^ 
results in a flux attenuation a factor of ~ 40 smaller at 1 pm 
than at 5100 A (in the rest-frame). The case is even more 
extreme for Compton-thick AGN with their expected hy- 
drogen column de nsities in excess of A^h > lO'^* cm~^ (e.g., 
iDaddi et al.ll2007h . Z. 

As we have shown in iLandt et al.1 (|201ll ). AGN black 
holes masses can be estimated from the near-IR virial prod- 
uct based on the 1 ^m continuum luminosity and the width 
of the Pa/3 (or Paa) broad emission line because the widths 
of the broad Paschen li nes are well-correla ted with those of 
the broad Balmer lines l|Landt et al.ll2003 ). The advantages 
of using the near-IR instead of the optical virial product 
are threefol d. First, since Paa and Pa.p are observed to be 
unblended (|Landt et al.ll200^ '). their width can be reliably 
measured. By contrast, the II/3 broad emission line is gen- 
erally observed to be strongly blended with both Fe II and 
He II A4686 and often shows a "red shelf" most likely formed 
by weak Fe II multiplets and He I. Secondly, the AGN con- 
tinuum around 1 pm is free from major contaminating com- 
ponents and may be easily determined, unlike the optical, 
which can suffer from an Fe II pseudo-continuum. Finally, 
since the near-IR is much less affected by dust extinction, it 
can potentially be applied to dust-obscured AGN. 



4 CONCLUSIONS 

We have presented the near-IR radius-luminosity relation- 
ship. Using near-IR spectra of 14 reverberation-mapped 
AGN obtained through a slit small enough that it excludes 
most of the host galaxy starlight, we have fit the relation- 
ship between the radius of the H/3 broad-emission line region 
and the integrated 1 /im continuum luminosity. The best-fit 
slope is in most cases ~ 0.5 ± 0.1 and always consistent 



with a value of 0.5, which is expected based on simple pho- 
toionisation arguments. The near-IR R-L relationship as an 
alternative to the optical relationship is expected to be rel- 
evant particularly for dust-obscured AGN. Black hole mass 
estimates can be obtained from the near-IR virial product 
based on the integrated 1 fim continuum luminosity and the 
widths of the strong and unblended Paschen broad emission 
lines Paa or Pa/3. 
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